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Abstract 

We present analytical expressions for the time-dependent and stationary 
probability distributions corresponding to a stochastically perturbed one- 
dimensional flow with critical points, in two physically relevant situations: 
delayed evolution, in which the flow alternates with a quiescent state in which 
the variate remains frozen at its current value for random intervals of time; 
and interrupted evolution, in which the variate is also re-set in the quiescent 
state to a random value drawn from a fixed distribution. In the former case, 
the effect of the delay upon the first passage time statistics is analyzed. In 
the latter case, the conditions under which an extended stationary distribu- 
tion can exist as a consequence of the competition between an attractor in 
the flow and the random re-setting are examined. We elucidate the role of 
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the normalization condition in eliminating the singularities arising from the 
unstable critical points of the flow, and present a number of representative 
examples. A simple formula is obtained for the stationary distribution and 
interpreted physically. A similar interpretation is also given for the known 
formula for the stationary distribution in a full-fledged dichotomous flow. 

Keywords: Randomly interrupted flow, critical points, time-dependent 
distribution, stationary distribution. 
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1 Introduction 

Physical systems are invariably subject to stochastic perturbations of various 
kinds. The effects of noise upon the evolution of classical dynamical systems 
have been studied quite extensively for several decades now ||. While 
much of the work in the area has been concerned with Gaussian noise (which 
offers several analytical advantages), other forms of noise are also relevant 
in specific instances. The inclusion of a finite correlation time for the noise 
(a closer approximation to physical reality) has also led to the uncovering 
of several interesting features that are absent in the idealized case of delta- 
correlated noise. 

In many physical situations, the stochasticity affects the dynamics more 
strongly than as just a small perturbation of the deterministic evolution. An 
important case is that of random interruption of the deterministic flow of a 
system variable x; during each interruption (the 'quiescent' state), x either 
remains frozen at its current value, or is re-set to a random value drawn 
from a prescribed distribution, and remains frozen there. Further evolu- 
tion (the 'active state') is resumed after a random interval of time. We call 
these possibilities delayed evolution and interrupted evolution, respectively. 
A main purpose of this paper is to obtain exact expressions for both the time- 
dependent and the stationary probability distributions of the driven variable 
x in both these situations. While problems of this general type have been 
studied right from the early days in a variety of specific contexts ranging 
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from chromatography || to the collision broadening of line shapes § and 
other relaxation phenomena 0, our focus here is on an altogether different 
aspect: the elucidation of the role of the stable and unstable critical points 
of the flow in determining these distributions. 

The random interruptions are most simply modeled by assuming that 
they are switched on and off according to a Poisson pulse process of intensity 
A 0. The mean duration A -1 of the active and quiescent states provides 
a time scale whose interplay with the characteristic time scale of the deter- 
ministic evolution will be of interest. (For simplicity we assume that the 
switching between the two states occurs at a common mean rate A. All our 
results can be easily extended to the case of unequal switching rates in a 
straightforward manner.) We identify the active and quiescent states by the 
state labels £ = +1 and £ = — 1 respectively. It is useful to note that the 
random variable is in fact a symmetric, stationary dichotomous Markov 
process (DMP) with a correlation time r c = (2A) _1 . The conditional prob- 
ability density P(x, £, t\x , £ ) will occasionally be abbreviated to P±(x,t) 
when it is convenient to do so without causing confusion. 



2 Delayed evolution 
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2.1 Solution for P(x, t\xo, £o) 

Let the deterministic evolution (in the active or £ = +1 state) be given by 

± = /+(*)■ (!) 

The first step is to identify the critical points of the flow given by Eq. ([!]), 
i.e., the roots of f+(x) = 0, and the direction in which the flow occurs in each 
interval between successive critical points. Clearly, for an initial condition 
x(0) = Xq lying between two such points, x(t) is restricted to this interval for 
all t > 0. We therefore expect P±(x,t) to involve step functions that reflect 
this fact. 

The device we use to solve the problem of obtaining the exact time- 
dependent probability distributions for delayed evolution is to map the prob- 
lem onto that of a particular dichotomous flow. Since x vanishes identically 
in the state £ = — 1, the evolution equation for x at any time t can be written 

as 

i = (I±M) /+w . (2) 

Equation (g) implies that {#(£),£(£)} constitute a stationary Markov process. 
Denoting the anti-derivative of l/f + (x) by q(x), i.e., q(x) = J dx/f + (x), we 
see that q is a DMP that switches between the values 1 and 0. Therefore the 
(non-stationary!) process y = q — |t is given by the stochastic differential 
equation y = ((t) where ((t) is a DMP switching between the values \ and 
— ~: in other words, y is a persistent diffusion process []TTf| . Its conditional 



densities P±(y,t) (retaining the symbol P for these, in a slight abuse of 

notation) therefore satisfy the telegrapher's equation 
( d 2 d 1 d 2 \ 

The general solution of Eq. (|3|) can be written down |L2| (using, for in- 
stance, the method of characteristics), once the initial conditions P±{y,0) 
and P±(y,0) are specified. To obtain the latter, we begin with the initial 
conditions on the original densities, namely, 

P(x, ±, 0|z„, ±) = S(x - x ) , P(x, ±, 0|x , T) = • (4) 

These translate into the conditions 

P(y, ±, 0\y , ±) = 5(y - y ) , P(y, ±, 0\y Q , =f) = (5) 

for the densities of the y-process, where yo = q(xo). The initial conditions 
on the time derivatives of the densities are found by recalling that P + (y,t) 
and P_(y,i) actually satisfy the coupled first order equations (d t ± ~d y + 
\)P±(y,t) = \P T (y,t). Insertion of Eqs. (|) in the latter leads to the initial 
conditions 

P(y,±,0\y ,±) = ^S'(y-y )-XS(y-yo) , P(y, ±, 0\y , t) = A%-y ) • 

(6) 

Using Eqs.(^|) and (||) in the general solution of Eq. (Q) and simplifying, we 
finally arrive at the solutions given in Eqs. (|T0|)-(|1^) below for the set of 
densities P(x, t\x , £ ). Let 



T(x ,x) = / dx'/f+(x') = q(x) - q(x ) (7) 

JXQ 
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denote the time taken to move from xq to a given point x under the (au- 
tonomous) flow of Eq. ([]]). Correspondingly, the location of the point reached 
under this flow in a given time interval t, starting from xq, is given by the 
expression 

q~ 1 (q(x ) + t) = x(x Q ,t) . (8) 

In reverting from y to the original variable x, use is made of the relationship 

8 (t — T(x , x)) 

H*-xM)= ]f+ l; " ■ («) 

Then, omitting the arguments in X(xo,t) and T(xq,x) for brevity, we find: 



P(x,+,t\x ,+) = e 



-At 



S(x-X) + 



A / T \ 2 



+ jj^ [j—f) J i ( 2A y/T(t - T) ) 9(x - x ) 9(X - ^0) 



P(x,+,t\x ,-) = P(x,-,t\x ,+ 

Xe- Xt 



\U(x)\ 

and 



J 2XJT(t - T) 9{x - x Q ) 6(X - x) ,(11) 



-XL 



5(x - Xq) + 



A /T — A 5 



P(x,-,t\x ,-) = e 

+ jj^ {^) h [2Xy/T(t - T) j 6(x - xo) 9(X - ^2) 

where I n is the modified Bessel function of order n. The step functions in Eqs. 
(|l0l)-(|T2"|) refer to the situation in which f+(x) > in the interval concerned. 
They originally appear as 6{T) 6{t — T), so that when f+(x) < the product 
of step functions in these equations is replaced by 6{xq — x) 9(x — X). 



2.2 First passage time distributions 

An interesting question that arises naturally in the case of evolution that is 
stochastically delayed as above is the following: what is the effect of the delay 
on the corresponding level crossing (or first passage time) statistics? As the 
flow does not reverse direction in any interval between two successive critical 
points, it is clear that x(t) is a non- decreasing [respectively, non-increasing] 
function of t for all t > if f+(x) > [respectively, f+(x) < 0] in the 
interval concerned. It follows at once that P(x, +, t\xo, £o) dx is also equal to 
the probability Q(t, x\xo, £o) dt of reaching the point x for the first time at 
time t. Thus the quantity |/+(x)| P(x, +,t\xo,£o) = Q(t,x\xo,£o), expressed 
in terms of the appropriate variable (namely, t), is the first passage time 
density for reaching the point x, starting at t — from the point xq and in 
the state £o- With the help of the useful relationship in Eq. (£5p we obtain, 
after a bit of simplification, 



Q(t,x\x , +) = e 



-At 



d 



5(t-T) + -I [2Xy/T(t-T))9(t-T) 



(13) 



and 



Q(t, x\x 0} -) = Ae- At / 2\JT(t - T) 6(t - T) . (14) 



The corresponding generating functions (i.e., Laplace transforms) are given 
by relatively simple expressions: 

Q(s, x\x , +) = e-^ T e A2T /( s+A ) (15) 
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and 

Q(s, X \x ,-) = I ^e^ T e^^. (16) 

As these expressions reduce to unity when s — 0, it is evident that the first 
passage time distributions concerned are properly normalized. 
The mean first passage time is then 

(t{x -> x))+ = 2 T(x , x) , (t(x -> x)>_ = 2 T(x , x) + A^ 1 , (17) 

where the subscripts on the angular brackets indicate the starting state. 
The extra A -1 in (t)_ is understood as follows. We expect to exceed 
(t) + by just the mean forward recurrence time in the £ = — 1 state; and 
since the DMP is governed by an uncorrelated Poisson process (implying an 
exponential waiting time density), this is the same as the mean duration of 
this state, i.e., A -1 . 

The mean squares are 

(t 2 (x -> x))+ = 4T 2 + 2TA- 1 , (t 2 (x = 4T 2 + 6TA" 1 + 2 A~ 2 . 

(18) 

The stationary a priori probabilities of being in the active and quiescent 
states being equal, the net mean first passage time and its mean square are 
therefore 



(t(x -> x)) = -({t)+ + = 2T( l0l i) + (2A)~ X =2T(x ,i) + r c (19) 
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and 

(t 2 (x -> x)) = \{{t 2 ) + + (t 2 >_) = (2 T(x , x) + A" 1 ) 2 = 4 (T(x , s) + r c ) 2 , 

(20) 

where we recall that r c = (2A)" 1 is the correlation time of the DMP 
controlling the switching between the evolving and quiescent states. The rel- 
ative fluctuation in the first passage time (the ratio of its standard deviation 
to its mean) is therefore given by 

Jt c UT + 3t c ) 
AtOro - s) = L . (21) 

I 1 + T c 

We note that At tends to a constant value (= ^3) as T(x ,x)/t c — > 0, while 
it decays like (T(x , a;)/r c ) _ 5 for very large values of this ratio. 

These results can easily be generalized to the case of unequal mean du- 
rations of the evolving and quiescent states. 

3 Interrupted evolution 
3.1 Solution for P(x, £, t\xo, ^o) 

We turn now to the case when x is re-set instantaneously to a random value 
drawn from a prescribed normalized distribution <f>-(x) whenever a transition 
occurs from the active or evolving state to the quiescent state, and remains 
fixed at that value till the active state occurs again. It is evident that the 
re-setting can be taken to occur at any instant during the quiescent state, 
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without affecting the results. We note that x may be re-set randomly any- 
where in its total range (—00, 00), and not just in the interval between the 
two successive zeroes of f+(x) that contains its initial value xq. 

Once again, it is helpful to regard the switching between states as being 
triggered by the DMP £(t). Denoting the transition probabilities for this 
DMP by p(C,t\Co), we have 

p(±,*|±) = e~ A *coshAt , p(±,t|=F) = e _A 'sinhAt . (22) 

Each of these probabilities tends to \ as t — > 00. We therefore have the limit 



lim P(x, -, t\x , f ) = \ <P-{x) , (23) 

i— >oo Z 



where 

/oo 
<P_(x)dx = 1. (24) 
-00 

Similarly we may ask whether 



lim P(x, +, t\x , f ) = 7: <f>+(x) , (25) 



where (f>+(x) represents the stationary probability density of x in the active 
state. The conditions under which a normalizable, extended density <f>+(x) 
exists will be examined in detail in the next subsection. 

Owing to the instantaneous re-setting of x in each occurrence of the qui- 
escent state, it is no longer useful to map the problem to a stochastic differ- 
ential equation. However, it is clear that it is only the last such re-setting 
that controls the location of the final point x at time t, in each realization of 
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the random process. The exact expressions for P(x, £, t\xo, £o) are therefore 
easily found by enumerating the processes contributing to the probabilities 
and summing over the corresponding propagators. As £(£) is governed by a 
Poisson process of intensity A, the probability of a flip in occurring in 
an infinitesimal interval 8t is X5t, while the probability of the persistence of 
either state of £ for an interval t is exp(— At). Using these facts, and starting 
with the simplest case, we find 

P(x,-,t\x ,+) = E / ^y^^A^^^'^-W 

n=0 J0 l Zn J- 

= e~ A 'sinhAt0_(a:) . (26) 

While this is of course identical to <f>-(x) p(—, t\+) (and has no xo-dependence) , 
as one would expect, it is interesting to note that P(x, — , t\xo, — ) is not sim- 
ply equal to (f>_(x)p(—,t\—). Owing to the zero-transition contribution we 
obtain, instead, 

P(x,-,t\x ,-) =e- xt [5(x-x ) + (f)_(x)(cosh\t-l)] . (27) 

The non-trivial cases (corresponding to a final state = +1) can be ana- 
lyzed similarly, to obtain the following results: 



P(x,+,t\x ,+) = e 



-At 



S(x-X(x ,t)) 



roo ft 

' , (28) 



/oo rt 
dx' / dt' <p„{x') 5(x - X(x',t- t')) sinhAt' 
-oo JO 



and 



/oo rt 
dx' / dt' <p.(x) 5(x-X(x',t-t')) cosh At' .(29) 
-oo JO 
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Here X(x',t — t') = q~ 1 (q(x') + 1 — t'), in accordance with the definition in 
Eq. (D). The xo-dependence of P(x, +, t\xo, — ) is implicit in Eq. (|29"D ; it 
arises from the 'one-transition' contribution to P(x, +,t\xo, — ), and is given 
by A t e~ xt 0_(s o ). It is easily verified that the expressions in Eqs. (p6|)-(p9|) 
are correctly normalized, according to 

/oo /*oo 
P(x,±,t\xo,±) dx = p(±,t\±) , / P(x, ±,t\xo,^f) dx = p(±,t\^f) . 
-oo J — oo 

(30) 

The integration over x' in Eqs. (^) and (|29|) is obviously constrained by the 
5-function in each integrand: for a given x, it is restricted to those points 
lying on the trajectory X(x' ,t — t') that are carried by the flow of Eq. ([!]), in 
a time ranging from to t, into the point x at precisely the instant t. Using 
the ^-functions to carry out the integration over t' rather than x' helps us 
express the results in a suggestive form. Equation ([]) implies that 

6 {x -X {x \ t -t>)) = S ( t '-^(*'.*)) . (31 ) 



Using this in Eqs. (|^) and (|29|) , we finally get 

S(x-X(x ,t)) 

A 

+ 



P(x,+,t\x ,+) = e~ M 

A 



/ dx' (f)-(x') sinh A(t — T(x', x)) 



(32) 



/ + (ar) yx(a!,-t) 
and 

P(x,+,t|x ,-) = f X dx'(j)(x') coshA(t-T(x',x)) . (33) 

f+(x) Jx(x,-t) 

We note that it is /+, rather than |/ + |, that appears in these final expres- 
sions: The limits of integration as they stand here automatically ensure the 
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positivity of the probability densities concerned, because X(x, —t) lies to the 
left [right] of x if /+ is positive [negative] in the relevant interval, for t > 0. 



3.2 Stationary distribution in the active state 



It is trivially seen from Eqs. (p6|) and fl27|) that lim^oo P(x, — , t\xo, £o 



^ </>_ (x) as required. The interesting question is whether there exists a non- 
degenerate stationary distribution <f>+(x) in the active state, in accordance 
with Eq. (p5j ). In other words, what is the eventual outcome of the competi- 
tion between the attracting fixed point (s) towards which x evolves each time 
it is in the active state, and the random re-setting it undergoes whenever it 
falls into the quiescent state? 

The more familiar approach to this question would be to start with the 
differential equation for P(x, +, t\xo, £o), namely, (d t + d x /+ + A) P+(x, t) = 
A P-{x, t), in the limit t — > oo: the stationary distribution (/>+(x) must satisfy 
the first order equation 

\^f+(x) + Xj<f> + {x) = <f>.(x). (34) 

Moreover, it must be non-negative and normalizable. Equation (|34]) does 
make it clear that </>-(x) acts as the 'source' for <fi + (x). However, the dis- 
advantage of this approach is that a constant of integration is explicitly 
involved: if <j)+(x) = 4> + at x = x, then the formal solution of Eq. Q3~4j) reads 



4>+{x) = + e-r^' + [ X dx' ( ^441 e-CH*"W (35) 

Jx \f+{ X )J 
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where F(x) = (f+(x) + X)/ f+(x). It is not immediately clear how to choose 
x and determine 4> + , and how the explicit dependence on these quantities 
disappears in the final result for + ( X^j j clS it must. Presumably, the other 
requirements on (j)+(x) such as its normalizability play a role here. As this 
is an interesting (and instructive) question in its own right, we shall return 
to it shortly. Before doing so, however, we point out that the exact time- 
dependent solution in Eq. Q52D (or Eq. fl3"3])) already provides us with the 
following complete answer, circumventing these problems. 

Let a, j3 denote two successive zeroes of f+(x), and let the flow x = f+(x) 
be directed from a to j3 in the interval (a, (3). As a is a repellor, we have 
X(x, — oo) = a for x G (at, (3)- Therefore, passing to the limit t — > oo in Eq. 
(|D or (H), we obtain 

<P + (x) = -A- f tfe' 0_(x') e ~ xnx '' x) , s e (a, P) . (36) 

The same expression is valid in each such interval between successive critical 
points; the integration over x' runs from the repellor in that interval to the 
point x. Equation (|36| ) not only specifies the exact solution for the stationary 
distribution in the active state, but also helps determine when such a solution 
exists, whether it is normalizable or not, when a pile- up or a divergence of the 
density occurs at the attracting fixed point, and when there is a "leakage of 
probability" into the absorbing state. As all these aspects are best illustrated 
by considering a set of typical examples as case studies, we now turn to these. 
In what follows, 7 is a positive constant denoting a characteristic time scale 
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of the deterministic dynamics given by x = f+(x). 

(i) Stable drift: Consider the simple case f+(x) = — jx, corresponding to a 
simple attractor at x — 0. For definiteness, let us confine our attention to 
the region x > 0. The repellor in this region is x = +00, so that the formula 
of Eq. ( p6[) yields in this case the exact solution 

M*) = ^—] g J^dx . (37) 



On the other hand, the formal solution given by Eq. fl3"o] ) yields 



x- 1+x h 7 y s (x') A/7 

We require that <f)+(x) be normalizable (i.e., integrable). However, the factor 
x _1+A /7 leads to a divergence at +00 (the repellor), unless the factor in 
parentheses in Eq. ( [38|) vanishes as x — > +00. Therefore we must impose the 
condition 



x- 1+x h 7 (x') A/7 
which not only removes the divergence, but also eliminates all dependence 
on x and 0+ and leads to precisely the normalizable density given by Eq. 



(ii) Unstable drift; power law tail : Now let f+(x) = +7X, corresponding 
to an unstable critical point at the origin. We consider the region x > as 
before. As x = is now the repellor, Eq. (RBI) yields 



4>+{x) = — * /Y(z') (x') xh dx> . (40) 



7a ,l+A/ 7 J Q 
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On the other hand, applying Eq. d35|) we obtain 

<j> + {x) = -j^j- U+ x 1+A/7 + - f <t>-{x') (*') V7 dx') . (41) 

The problem of normalizability now arises from the factor x~ 1 ~^ 1 , and again 
occurs at the repelling critical point (here, at x — 0). To eliminate it, the 
factor in parentheses in Eq. ( f4"T|) must be required to vanish as x — > 0. The 
condition to be imposed is therefore 

+ x 1+A/7 = - f <j>_{af) (x') A/7 t&' , (42) 
7 Jo 

leading once again to the correct solution as given by Eq. ([10). 

The following point is noteworthy. An inspection of Eq. fliOD shows 
clearly that, in this case, the density <fi+(x) always has a power law tail. 
(When 4>-{x) is an even function of x, replacing x with \x\ in Eq. (|40|) yields 
(j)+(x) in the entire range — oo < x < oo.) The random re-setting of x in the 
quiescent state does compensate for the unstable drift towards infinity in the 
evolving state, but only to the extent of producing a normalizable stationary 
distribution in the latter state. It is insufficient to prevent a slow (power-law) 
fall-off of the latter, no matter how rapidly the fixed density 0-(x) falls off 
for large \x\, or even if its support is compact. Even then, the variance of 
x diverges unless A/7 > 2, i.e., the switching rate is greater than twice the 
drift coefficient. 

This example also represents a very direct way in which simple dynamics 
can generate the whole family of stable (Levy) distributions. As the (cumula- 
tive) distribution function of x has a tail ~ x~ x ^, a (suitably re-scaled and 
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shifted) sum of such variates, independently distributed, will have a stable 
distribution with exponent A/7, going over into a normal distribution for 



(hi) Higher order critical point; 'leakage' of probability : Next, consider the 
case f+(x) = 7X 2 . The flow is towards x = for all negative x, while it is 
directed away from this point and towards +00 for positive x. As before, 
the exact solution for (f>+(x) may be found by using Eq. (|3lf), or by using 
Eq. (^) and imposing the conditions required to eliminate divergences and 
ensure normalizability (at — 00 for negative x, and at x = for positive x). 
The result is, for x > 0, 



For x < 0, the lower limit of integration is —00 instead of 0. 

As in the preceding case, <fi+{x) has a power law tail. But the degenerate 
critical point at the origin has an even stronger effect. Any positive initial 
value Xq reaches +00 in a finite time (7x0) _1 under the flow in the evolving 
state. We may therefore expect some sort of "absorption" at +00: in terms 
of probability distributions, this would show up as a "leakage of probability" 
leading to a deficit in the total probability. This is indeed borne out: we find 
that 



A/7 > 2. 




(43) 




(44) 



but 




(45) 



18 



This loss in probability measure can also be given a simple interpretation. 
The probability of being injected into (x, x+dx) under the random re-setting 
is <p-(x)dx. The time taken to reach +00 from any x > under the flow is 
(7a;) -1 . The probability of remaining in the evolving state for this duration is 
exp(— A/72;). Therefore the total loss in measure is 4>-(x) exp(—\/ r yx)dx. 
In general, such a loss of measure occurs whenever the absorbing state at the 
attracting critical point is reached in a finite time under the deterministic 
flow, as will become clear from the examples to follow. 

It is also an instructive exercise to "unfold" the degenerate critical point 
at x = by starting, for instance, with f + (x) = 7x(x— e), and then examining 
how various quantities behave as e — > 0. 

(iv) Periodic boundary conditions; single critical point : In many physical 
applications, x is an angular variable, so that its range is compact, and 
periodic boundary conditions apply. Consider the case x = 7 | sinx|, with 
a fundamental period equal to n. There is only one critical point, and the 
flow is directed from towards tt. Given that <p-(x) = <ft^(x + tt), we seek a 
solution (j)+(x) that has the same periodicity property. Using Eq. (|36|), we 
find 

, / n A r x ft&nx'/2\ A/7 , , , . p . , . 

( f )+ ( x ) = ^— '— 0_ ( x > )dx' , xe 0,tt , 46 

7Siniio \ta.nx/2J 

together with (j)+(x) = <p+(x + it). The time taken to reach 7r from any x in 
< x < tt is infinite (although the range is finite, the velocity f+(x) vanishes 
as x — > it). Therefore there is no loss of measure in this case. However, a 
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pile-up of the density can occur at x = it, depending on the value of the 
ratio A/7: for a uniform </>_(x), for instance, we can show from Eq. ( [46|) that 
<f)+{x) is finite at x = n only if the switching rate A > 7. As x j 7T, we find 
(f>+(x) ~ — In (n — x) for A = 7, and 4>+(x) ~ (7T — x) _1+A / 7 for A < 7. Thus 
is divergent (though integrable) at the attractor, unless the switching 
rate exceeds the drift coefficient. 

(v) Periodic boundary conditions; degenerate critical point : As the periodic 
analog of Case (iii) above, consider f + (x) = 7 sin 2 x with the fundamental 
interval [0,7r). We find the periodic solution 

6+(x) = — ^ e (A/7)cota: r cj)Jx')e' (xh)cotx ' dx' , x 6 [0, it) . (47) 
7 sin x Jo 

Again, as the time taken to reach 7r from any x in < x < tt is infinite, 
there is no loss of measure owing to absorption at n. Moreover, now there is 
no divergence of <f>+(x) at x = it, either, because the velocity f+(x) vanishes 
more rapidly than linearly as x — » n. 

A flow like x = jx with periodic boundary conditions ([0, 1] being the 
fundamental interval, for instance), will produce a leakage of probability into 
the (finite) absorbing point at x = 1: now the velocity does not vanish as 
x | 1. It is easily shown that f$ <f> + (x) dx = f$ <f>-{x) (1 — x x ^) dx in this 
case. 

(vi) Periodic boundary conditions; two critical points : Finally, let us consider 
the case x — 7 sinx, with a period equal to 2n. The velocity f+(x) changes 
sign in the fundamental interval, and there is a repellor (at 0) as well as 
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an attractor (at 7r) in it. In the range < x < it, the solution for (j>+(x) 
is exactly as in Eq. fl46|) ; for tc < x < 2tt, the lower limit of integration is 
replaced by 2tt. The divergent behavior found as x | tt when A < 7 in Case 
(iv) above now occurs symmetrically on both sides of x = n, as one would 
expect. 



4 Stationary distribution in dichotomous flow: 
some remarks 



The form of the stationary distribution in Eq. (j36|) and the discussion fol- 
lowing it suggest a direct physical interpretation (amounting, in fact, to a 
simple heuristic derivation) of the well-known formula [13| , ]14[ for the sta- 



tionary distribution that obtains (under certain conditions) in the case of the 
full-fledged dichotomous flow given by the stochastic differential equation 

x = f{x) + g(x)£(t) . (48) 

Such a flow corresponds to the random alternation of two deterministic flows, 
given by x = f+(x) and x = f-(x), respectively, where f±(x) = f(x) ± g(x). 
It serves as a model for numerous physical phenomena (see, e.g., Refs. Hl4|] , 

The time-dependent probability densities P(x, £, t\x , £ ) (denoted by P±(x, t) 
for short) now satisfy the coupled first order equations (d t +d x f±+X) P±(x, t) = 
XP T (x,t) . In general, it is not possible to obtain a partial differential equa- 
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tion of finite order for the total probability density P(x,t) = P + (x,t) + 



P-(x,t): this possibility is restricted to certain special cases |ll|-[[18]. How- 
ever, if the deterministic flow is "dynamically stable", the stationary density 
P st (x) = lim^oo P(x, t) is found to satisfy the first order ordinary differential 
equation 

({f + g)(l- f -)P st ^ - 2\ f -P st = . (49) 

Here, the term "dynamic stability" denotes precisely the situation of interest 
in the present context: the existence of two different stable critical points in 
the flows x = f-\-(x) and x = f-(x) respectively, with no other critical point 
of either flow in between them. What is sought is the stationary distribution 
to which x settles down under the competition between the two attractors. 
(Thus the second attractor takes on the role played in the preceding section by 
the random re-setting of x.) A typical case that helps visualize the situation 
is as follows: a repellor at x = j3 and an attractor at x = b in the flow f+{x); 
and an attractor at x = a and a repellor at x = a in the flow f-{x). If 
(3 < a < b < a, the interval [a, b] acts as a trapping region supporting a 
non-trivial P st (x). The latter is given |13|,|14]| by the solution of Eq. (|49|), 
namely, 

P-(*)=A/--^exp( 2 A/^) , (50) 

where M is the normalization constant. We give a simple interpretation of 
this formula along the lines of that presented for Eq. (|36|) above, which 
suggests how the expression in Eq. ([5(]) can virtually be written down by 
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inspection. 

Consider, first, what happens if the system alternates regularly between 
the two flows, with a constant duration r in each state. Let us denote the 
solution trajectories in the two flows by X + (x,t) and X_(x,t), respectively. 
Thus, if we start in the + state from the point xo, we have x(t) = X + (x , r), 
x(2t) = X-(x(t), r), and so on. Setting x(2nr) = r\ n and x((2n + l)r) = ( n 
(where n — 0, 1, . . .), the dynamics is given by the two-dimensional map 



In the "dynamically stable" situation referred to above, there exist stable 
fixed points r]* and C* such that 77* = X_(X + (rj*)), (* = X + (X_((*)), im- 
plying a stable period-two cycle into which the original variable x(nr) falls 
asymptotically as n — > oo. Correspondingly, in continuous time the trajec- 
tory x(t) zig-zags between the values r/* and (*, alternately following the two 
evolution rules. An invariant density for x can therefore be defined in each 
of the two distinct alternating segments. Calling these (f>+(x) and (f>-(x), 
respectively, it is evident that (f>±(x)\dx\± = r^dt, or <f>±(x) = (r |/±(x)|) _1 , 
i.e., essentially just the reciprocal of the corresponding Jacobian. As the + 
and — states are equally probable, the invariant density of x becomes 



Taking into account the fact that the flows are in opposite directions, i.e., 
that f+(x) and f-(x) have opposite signs, the quantity in parentheses in Eq. 



(51) 




(52) 
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(|52|) is proportional to fz l — /+ 1 = g(g 2 — / 2 ) -1 . This explains the origin of 
the factor preceding the exponential in Eq. (|50|) . 

It remains to understand the exponential factor in Eq. (|50|). Returning 
to the case of random switching between the two flows, it is evident that the 
weight factor r _1 in Eq. ([52]), applicable when the switching occurs at regular 
intervals, must be replaced by a product of two probability factors, one for 
each of the two flows. This factor is the same as that already discussed 
in relation with Eq. (|36D in the case of interrupted evolution, namely, an 
exponential of the time required to reach the point x from the corresponding 
repellor - or rather, from the attractor in the other flow, recalling that the 
trapping region (the support of P st (x)) is bounded by these points. Thus, 
in the example mentioned above (in which (3 < a < b < a)), we have the 
factor exp[— AT + (a, x)] exp[— ATL(&, x)], where we have used subscripts on 
T to indicate the corresponding flow. Since T±(u,x) = dx' / f±(x') and 
/± = / ± 9, we find 



This yields precisely the exponential factor occurring in the solution given 
by Eq. (|50|) . The "mechanism" underlying this formula is therefore made 
more manifest by re-writing it in the form 



r. 




(53) 



p st (x) = J Y (\u(x)\- i + \f-(^r 1 ) e 



AT+(o,«) e -AT_(i),i) 



(54) 
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for the configuration of attractors under consideration, with obvious minor 
modifications for the other possible configurations. 



A simple example is provided by the exactly solvable case |L7| of linear 
dichotomous flow, x = -71 + c£(t), where 7 and c are positive constants. 
Here (3 = —00, a = — c/7, b = c/7, a = 00. The stationary density of x has 
support in the interval [—c/7 , c/7], and is proportional to (c 2 — 7 2 x 2 )~ 1+A / 7 . 
We can now readily understand its interesting features as follows. Since 
T + (x , c/7) and T„(x , —c/7) diverge, there is no loss of probability measure 
owing to the absorption at the attractors in this instance. However, as we 
may now anticipate (based on the remarks made in Case (iv) above), P st (x) 
is indeed divergent at the (finite) endpoints ±0/7 when the switching rate A 
is smaller than the drift coefficient 7. 
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